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Abstract 

O" 

CN Based on the external perturbation that disturbs the system only slightly from its equilibrium 

Oh- 

q_) ■ position we make the Taylor expansion of the pressure of a quark gas. It turns out that the first 



o 
o 



X 



term was used in the literature to construct a Hard Thermal Loop perturbation theory (HTLpt) 
within the variation principle of the lowest order of the thermal mass parameter. Various thermo- 



r^ • dynamic quantities within the 1-loop HTLpt encountered overcounting of the leading order (LO) 

q. contribution and also required a separation scale for soft and hard momenta. Using same varia- 

<— j , tional principle we reconstruct the HTLpt at the first derivative level of the pressure that takes 



into account the effect of the variation of the external source through the conserved density fluc- 
tuation. This modification markedly improves those quantities in 1-loop HTLpt in a simple way 



> 

i — . 

(^ , instead of pushing the calculation to a considerably more complicated 2-loop HTLpt. Moreover 



the results also agree with those obtained in the 2-loop approximately self-consistent ^-derivable 
Hard Thermal Loop resummation. We also discuss how this formalism can be extended for the 
higher order contributions 

PACS numbers: 12.38.Cy, 12.38.Mh, ll.10.Wx 
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I. INTRODUCTION 



The HTL resummation developed by Braaten and Pisarski [1] has been used to calculate 
various thermodynamic quantities in the literature based on two methods. The 2-loop 
approximately <l>-derivable approach was developed by Blaizot et al |2|, which produces a 
correct LO and plasmon effects for thermodynamic quantities (e.g, entropy density, number 
density etc.) and also for quark number susceptibility (QNS) J3j. On the other hand the 
HTLpt using variational principle through the lowest order of the thermal mass parameter 
was developed for pressure at the l-loop_level by Andersen et al 4J, which is at present 



151 



pushed to the 2-loop |5| and 3- loop level [6j . However, the 1-loop HTLpt pressure has a bad 
perturbative LO content, in the sense of severe over-inclusion of the effect of order g 2 as the 
HTL action is accurate only for soft momenta and for hard ones only in the vicinity of light 
cone. Such problem is cured (or at least pushed to higher orders) only after going to 2-loop 
level in HTlpt 5|, which is indeed a considerably more involved calculation. A very recent 
1-loop QNS calculation |7| from pressure in HTLpt |4J had, obviously, the problem of over- 
inclusion of the order g 2 . Moreover, it required an ad hoc separation scale to distinguish 
between hard and soft momenta. On the other hand the calculation of QNS in Ref. [8] 
dealt with the imaginary part of the charge-charge correlator in the vector channel and 
required to show the charge conservation. It also encountered some technical difficulties and 
over-inclusion problem in order g 2 as discussed in Ref. J3J. Also the Landau damping (LD) 
contribution was discussed but ignored. 

The equation of state (EOS) of strongly interacting matter at nonzero baryon density 
and high temperature is a subject of great interest for wide spectrum of physicists. Also 
QNS is a topical quantity in view of the ongoing efforts towards understanding the actual 



nature of the QGP Q as QNS plays an important ro.e U U * locating the critical end 

point in QCD phase diagram. As it stands the LO thermodynamic quantities |2|, |4j and 
QNS [3|, LD, IsJ] in HTL approximation |l| led to different results. This requires a detailed 
analysis of the leading order quantities within the HTLpt before extending it to the higher 
orders. In view of this we do not aim at higher orders calculations, rather we intend in 
this article to sort out the problems in 1-loop HTLpt, which produced different results from 
that of the ^-derivable approach within the HTL resummation |l| , and finally arrive at a 
consistent result despite the use of different approaches. 



The paper is organised as follows. In Sec. II we briefly discuss some generalities on 
fluctuations, correlation functions and susceptibilities based on external disturbance to a 
physical system. In Sec. Ill we modify the HTLpt based on the external disturbance. In 
Sec. IV the HTL thermodynamics in presence of the quark chemical potential and then QNS 
in LO are obtained. We also checked, both numerically and analytically, the perturbative 
content of LO HTL QNS in Sec. V. Finally, we conclude in Sec. VI. 

II. GENERALITIES 

A. Fluctuation and Susceptibility: 



Let O a be a Heisenberg operator where a may be associated with a degree of freedom 
in the system. In a static and uniform external field T w the (induced) expectation value of 
the operator O a (0, x) is written 12[ as 



<o« (o. ^»^ = - ^f^luZ^ =vJ A <°« (o. *)> 



(i) 



where the translational invariance is assumed, V is the volume of the system and 1-L ex is 
given by 

n ex = -J2 [ d3x ° a (o. ~$) ?<* • ( 2 ) 

The (static) susceptibility Xaa is defined as the rate with which the expectation value 
changes in response to an infinitesimal change in external field, 



Xaa{T) 



dT„ 



T- 



= f d 3 x (o a (0,lt)O a (0,t) 
=o J ^ 



(3) 



where (O a (0,x)O a (0,0)) is the two point correlation function with operators evaluated at 
equal times. There is no broken symmetry as 



(0«(o,^))| 



^^o 



a(o,t; 



o 



T^0 



(4) 



B. Thermodynamic Relations: 



The pressure is defined as 



T 
P= v lnZ 



(5) 



where T is temperature, V is the volume and Z is the partition function of a quark- ant iquark 
gas. The entropy density is defined as 

dV 
S = £. (6) 

The number density for a given quark flavour can be written as 

_ dV 1 Tr [Afe-PW-^] (AT) 
P ~~d^~ V Tr[e-^«-M0] ~ ~V~ ' ^' 

with M is the quark number operator and \x is the chemical potential. If /i — > 0, the quark 

number density vanishes due to CP invariance. 

The QNS is a measure of the response of the quark number density with infinitesimal 

change in the quark chemical potential, p, + 5p, at p — Y 0. Under such a situation the 

variation of the external field, T a , in ([2]) can be identified as the quark chemical potential 

p and the operator O a as the temporal component (Jo) of the external vector current, 

J a (t,x) = ipT a ip, where r CT is in general a three point function. Then the QNS for a given 

quark flavour follows from (jHJ) as 

d 2 P 

li=0 



*m ^ I 



/ rf 4 x (j (0,f)J (0,0)) =-limRen^(0,p), 



M =0 ^ 

where the number operator, Af = J Jo(t, x) d 3 x = J ip(x)T ip(x)d 3 x and Uq (uj p ,p) is the re- 
tarded time-time component of the Fourier transformed vector correlator H av {oj p ,p) with 
an external momentum P = (u p , |p| = p). To write (jSJ) in such a compact form we 



have used the fluctuation-dissipation theorem and the quark number conservation 12|, [13] , 
lim^olmilo^u^p) oc 5{u p ). 



III. MODIFICATION ON HARD THERMAL LOOP PERTURBATION THEORY 



□ 



The HTL Lagrangian density for quark including HTL correction term [L4j is written as 

C-HTL = £-QCD + 5£hTL 

= ^i llx D^ + m\i) lvL ( -r^—p U , (9) 

where R = (1, r) is a light like four- vector, ip and ^ are the fermionic fields, D is the 
covariant derivative, () is the average over all possible directions over loop momenta. The 
second term is gauge invariant, nonlocal and can generate iV-point HTL functions [1], which 



are inter- related through Ward identities. Now m q , is the quark mass in a hot and dense 
medium, which depends on the strong coupling g, temperature T and the chemical potential 
ft. Despite these facts m q is treated in (^ as a parameter much like the rest mass of a quark 



H 



and a HTLpt has been developed |4( around this rest mass (i.e., m ) by reorganising the 
HTL terms where m 2 was treated as the order of (gT)° for a hot system. In this way the 
effect of m 2 is taken into account in higher orders much like a variational principle. For a hot 
and dense system we will also treat the mass parameter as the order of (gT)° and (g^) , and 
reorganise the HTL term based on the variation of the external source and setting it zero 
at the end. In this way the effect of m 2 is taken into account to the higher order variations 
of the external source and thus to the response of the system. 

We note that the covariant derivative usually contains background field or any source, 
j depending upon the physical requirement. To motivate the perspective we define the 
covariant derivative D^ as 

D" = [V - iS"°(j + Sj)] = [D» - i5»° Sj] . (10) 

We note that T> contains gauge coupling and _D M = T>^ — iS^j, and 5j is an infinitismal 
change in external source to which the response of the system can be calculated, as discussed 
in Sec. II. Later it can be identified with a variation of some physical quantity depending 
upon the requirement of the system under consideration. 
Now, expanding the second term in ([9]), we can write as 

ChtlU + Sj) = Jp ($ + E) ^ + Sj ijT ij + 5f i^f^j + O(df) 

= ChtlU) + 5 J #W + Sj 2 i^f^ + O(Sf) (11) 

where the various iV-point functions in coordinate space are generated as 

£ = ml (-*J) , r = *- 7 „ - ml (*%£-) , L 00 = 2^ (fi^\ ,(12) 
q \iR-D/ M q \(iR-D) 2 / q \ (iR-Df /' V ; 

where these functions can easily be transformed into momentum space [15| . We now note 
that these iV-point HTL functions in ffT2|) are also inter-related by Ward identities. In HTL- 
approximation the 2-point function, S ~ gT (quark-self energy) is of the same order as 
the tree level one, Sq 1 ^) ~ IjC ~ gT (in the weak coupling limit g « 1), if the external 
momenta are soft, i.e., of the order of gT. The 3-point function is given by gT u = g^y+STy), 



where 5T V is the HTL correction. The 4-point function, g 2 T U(T , is higher order and does not 



exist at the bare perturbation theory and only appears within the HTL approximation 
Now considering the HTL Lagrangian in (ITTj) . we can write the partition function 



m- 



15| as 

Z[/3;j + 6j] = I V[i)\D[i)\D[A]e^ dixCHTL ^^ j+5l \ (13) 

where (3 = 1/T, is the inverse of the temperature and A is a background gauge field. 
The pressure can be written as 

V[(3:j + Sj} = ^lnZ[(3:j + Sj}, (14) 

where the four-volume, V = (3V with V is the three-volume. 
Expanding V in Taylor series around Sj one can write 

V\j3:j + 6j] = V[f3;j] + 5j V'[f3;j + 5j}\ 5 ^ + 6 -^- V"[P;j + 5j]\ 5j _ + .(15) 

The first derivative of V w.r.t. j is related to the conserved density in (00) whereas the 
second derivative is related to the conserved density fluctuation in Q. The above expansion 
in (I15p is very important for a resummed perturbation theory. We now note that a HTLpt 
was developed in Ref. [4J by considering the first term in ( TT5|) with j = 0, which caused an 
over- inclusion of the LO pressure. This was cured by going into two-loop level in HTLpt J5j, 
which is of-course a very involved in nature. As we will see below this could easily be 
corrected if one constructs a HTLpt at the first derivative level of V in ( JTBl where the effect 
of the variation of external field is taken into account. 
Now V can be obtained as 

dV[(3;j + 5j] = ^—- j ' V$]V[il>]V[A] J ' d'x ^(x)T \J]^(x)e(^ dixC —^^) (16) 



dj 



where we have used (TT2|) . The full HTL quark propagator in presence of uniform j can be 
written as 

. n( . / V$]V[il>]V[A]il> a (x)M*') ^P (i I d 4 xC HTL (^ ^ j)) 
^ J ^ X >-- $V$]V[ij)V[A)^(iSd±xC HTL ^^3)) ' [ } 

We now note that this full HTL propagator, S\j], is indeed difficult to calculate and we 

n n 

would approximate it by 1-loop HTL resummed propagator [1|, [14|, S*[j] and also other 
HTL functions below. 
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Now using (ITT)) and performing the traces over the colour, flavour, Dirac and coordinate 
indices in f TT6|) one can write 



dV[/3;j + Sj] 



Oj 



Sj=0 



d 4 K 



tr{S*{j}(K)T {j](K,-K;0)} 



(18) 



where 'tr' indicates the trace over the colour, flavour and Dirac indices. 
Similarly, we obtain V" as 



&r\J3;j + 6j] 



df 



-N c N f T 



d 3 k 



x ^Tr[^[j](A') T [j](K,-K;0) S^[j](-K) T [j](K,-K;0) 
ho 



-S*\j}{K)r 0O \j](K,-K-0,i 



(19) 



where Nj is the number of massless flavours, N c is the number of colour and 'Tr' indicates 
the trace over only the Dirac matrices. We have also used an identity based on unitarity of 
S*\j] as 

dS ™ K) = -S*V\W ^^f^ ^\j\{K) = S*\j\{K) r \j](K,-K;0) S*[j](K) (20) 

Now, if we identify j as the quark chemical potential /i, and Sj as its change 5fi, then 
(fT9l) would represent the QNS as 



Xtf) 



dp 
dpi 



d 2 V 



M-s-0 



dp? 



-N c N f T 



At^O 



d 3 k 



xJ2 Tr [S*(K)T (K, -K; 0)S*(K)T (K, -K; 0) - S*(K)T 00 (K, -K; 0, 0)] , (21) 

k =(2n+l)niT 

where the temporal correlation functions at the external momentum P = (u p , \p\) = 0, is 
related to the thermodynamic derivatives 1 . The first term in the second line of (121 j) is a 
1-loop self-energy whereas the second term corresponds to a tadpole in HTLpt with effective 
iV-point HTL functions. 

Now, ( TTHj) represents the LO net quark number density in presence of uniform external 
field p as 

d A K 



&P 
P(M = ^ 



(27T) 



rtr[5>](AT)r M(A:,-A:;0)] 



1 As already discussed Ref. [8( dealt with the dchnition of QNS that involves the static limit of the imaginary 
part of the dynamical charge-charge correlator. If one uses the number conservation directly, viz., the 
imaginary part of the charge-charge correlator is proportional to S(uj p ), then it becomes equal to (|21[) in 
which charge conservation is in-built by construction. 



N c N f T 



d 3 k 
(2*f 



J2 Tr[S*(K)r (K,-K;0)} 



(22) 



k =(2n+l)iriT+n 



d 




FIG. 1: 1-loop Feynman diagram in HTLpt for quark number density, p q that originates with the 
variation of \i of the 1-loop HTLpt pressure. The dashed line represents the background field. The 
solid blobs are 1-loop resummed HTL TV-point functions. 



d 
Tfjl 





FIG. 2: Same as Fig[T]but for the lowest order bare perturbation theory. 




FIG. 3: The 7^ sign indicates that it is not the correct diagram in the right hand side (rhs) as the 
ft variation is not taken into account properly. The diagram in rhs actually corresponds to p q that 
was obtained in Ref. 



u 



When jjl — > 0, the net quark density in (122]) would vanish as there is no broken CP 
symmetry, which becomes consistent with fl4]). Also, (122]) constitutes a LO HTLpt in the 
first derivative level of V (see Fig. [1]) similar to the usual perturbation theory where the 
bare iV-point functions (see Fig. [2]) are automatically replaced by the 1-loop resummed HTL 
iV-point functions. This suggests that HTL resummation technique provides a consistent 
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perturbative expansion if one goes beyond the lowest order perturbation theory. In contrast 
Ref. [7J did not employ the variation of the external source as done in (fl5|) . which leads to 
Fig. E] with a bare vertex for the calculation of the net quark number density. This resulted 
in overcounting of the LO QNS. It also required an ad hoc separation scale to distinguish 
between soft and hard momenta. 

Below we briefly outline some of the essential quantities in HTL resummation |l| , which 
are required to compute (1221 . The resummed HTL propagator in 1-loop approximation for 
momentum K is given as 

70 - 7 • k 70 + 7 • k 



S*(K) 



2D + (k ,k) 2D4k ,k) 



with 



D±(k ,k) 



k ± k 



m„ 



1 / k 

2 l 1=F lfc 



'ln^±£±l 



k — k 



(23) 



(24) 



m„ 



g_ 

6 



rp2 



fr 



7P 



(25) 



where g 2 = 47ra s , a s is the strong coupling. Now, the zeros of D± describe [17[ the in- 
medium propagation or quasiparticle (QP) dispersion of a particle excitation with energy 
u + having chirality to helicity ratio +1, and of a mode called plasmino with energy U- having 
chirality to helicity ratio — 1. In addition, D± contains a discontinuous part corresponding 
to Landau Damping (LD) due to the presence of Logarithmic term in (j2"4"j) . Using these 
general properties of the quark propagator one can obtain the in-medium spectral function 
for quarks. 

The pole part of the spectral function can be written as 



u) 2 - k 2 
g±{u,k) = * S(u 



0J±) + 



oj- 



k' 



2m 2 q 



8{ui + ui 



T) » 



(26) 



as D + has poles at uj + and — w_ whereas those of £)_ are at w_ and — uj + . 

For /cq < k 2 , there is a discontinuity in In ^z^ as lny = In \y\ — in , which leads to the 
spectral function, f3±(u), k), corresponding to the discontinuity in D±(k , k) as 



1 1 

/9±(w, k) = — Disc — — — — 
n D±{k ,k) 



— Im 

7T D±(ko,k) 



e-J-0 



5(±|-i)e(^- 



u 



WT *-^(±1 



2fc U1 fe- 



)"+ »5f(l=Ffi 



2 ■ 



(27) 



The zero momentum limit of the 3-point HTL function can be obtained from the Ward 
identity [15|, [l6| as 





T°{K, -K; 0) = — (S*- l (K)) = ai + lf?-k 



where 



with 



dk 



a±b = -D'±(k ,k), 



D' °- 



2m\ 



k =f k kl-k 2 



(28) 

(29) 
(30) 



IV. THERMODYNAMICS AND QUARK NUMBER SUSCEPTIBILITY 

We first obtain the net quark density p(T,fi), which is then used to obtain various ther- 
modynamic quantities, viz., pressure, entropy density and QNS. 

A. Free case 



In free case the number density can be written from (1221) as 

d 3 k 






(31) 

k =(2n+l)TriT+fi 

where the 3-point function is T = 7o and the 2-point function is the free quark propagator 
for momentum K is given as 

7o - 7 • k 70 + 7 • k 



S f (K) 



with 



2d + (k ,k) 2d_(fc ,fc) : 



d± = —k ± k . 



Using ( 13 2 p in ( 13 ip and performing the trace over Dirac matrices, we get 



p f (T,fi) = 2N c N f Tj 



d 3 k 
(2^f 



E 



1 



ko — k k + k 



(32) 



(33) 



(34) 



For eva 



tion 16] as 



fc =(2n+l)7riT+/i 

uating the frequency sum in (151]) . we use the standard technique of contour integra- 



1 
2~Ki 



c 



1 



1 



k — k k + k 



— tanh f I dkn = x 

2 I 2 J ° 2 2m 



-2iri) y^ Residues . (35) 
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It can be seen that the first term of (1351) has a simple pole at ko = k whereas the second 
term has a pole at fco = —k. After calculating the residues of those two terms, the number 
density becomes 



P f (T,n) 



-N c N f 



2N c N f 



d 3 k 

d 3 k 
(2^)3 



tanh 



P(k - fi) 



tanh 



P(k + ii) 



n(k — ji) — n(k + /i)] , 



(36) 



where n(x) = l/(e l3x + 1), is the Fermi distribution function. 

Now, the pressure is obtained by integrating the first line of (136]) w.r.t. fi as 

j— 3 [(3k + In (1 + e-M-A) + In (1 + e ~^ k+ ^)] , 



(37) 



where the first term is the zero-point energy that generates a usual vacuum divergence 16 ]. 
The entropy density in free case can be written from pressure as 

d 3 k 



flpf 
S'(T jA i) = — = 2N c N f 



dT ~"°" J J (2tt) 3 
P(k-(i) 0(k + n) 



ln(l 



-0(k-n) 



) + In (1 + e-^+M) 



+ 



e /3(fc-M) + 1 e /3(fc+M) + 1 

The QNS is obtained as 



N c N f 



7vr 2 T 3 fi 2 T 



45 



X f (T) = ^[p f AT,») 



4N c N f (3 



m=o 



d 3 k 

(2nf 



n{k)(l-n{k)) = N f T 2 



(38) 



(39) 



B. HTLpt Case 



Using f !23p . ( 12 8 p in ( 12 2 p and then performing the trace over Dirac matrices, the quark 
number density in HTLpt becomes 



P H1L (T^) 



2N c N f T 



2N r N f T 



C iV/. 



d 3 k 

d 3 k 
(2^f 



E 



k =(2n+l)iriT+[j, 



DL 



D'_ 



E 



2m\ 



k -k k + k k 2 -k 2 \D + D_ 



(40) 



Apart from the various poles due to QPs in (140]) it has LD part as D±(ko, k) contain Loga- 
rithmic terms which generate discontinuity for k 2 < k 2 , as discussed earlier. Equation (14"U]) 
can be decomposed in individual contribution as 



,Q p 



LD, 



f (T,ri=p<"(T,ri + p w (T,ri. 



(41) 
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1. Quasiparticle part (QP) 



The pole part of the number density can be written as 



^" = Ww5i/c 



(2vr) 3 2m J c , 



1 1 

+ 



2m\ 



k -k k + k k 2 -k 2 \D 
x — tanh d/co 



D_ 



(42) 



In general residues for various poles in third and fourth terms in ( l4"2"j) can be obtained as 



Res 



2m 2 1 

k 2 -k 2 rT± 



■1 ; 



Res 



ko=U)±,—UJzi 



2m 2 1 
k 2 -k 2 Dl 



(43) 



fco=±fc 



where D ± (fco = ±Jfe) = ±^. 

1. First two terms in ( I42p give the same contribution as free case in (I36p . 



2. The third term has four simple poles at ko = u + , —u^,k, —k. After performing the 
contour integration the third term can be written as 



2m J c , kl-k 2 D + 2 2 2 



+ tanh h tanh tanh 



— tanh 
P{k + /J,) 



P{u+ - //) 



(44) 



3. The fourth term has four simple poles at fco = cu_, — u + , —k, k. After performing the 
contour integration the fourth term can be written as 



2m l 1 P^P(k°-») 



1 



2m J c , k 2 - k 2 D_ 2 



tanh 



dkn 



P 



tanh 



P(u- - //) 



+ tanh ^±±A + tanh ^tlA _ tanh W ' + '£ 



(45) 



Using (136|) . (144|) and (145]) in (1421) one can obtain the HTL quasiparticle contributions to 
the quark number density as 



P Q "(T,IM) 



-N c N f 



2N c N f 



d 3 k 
J2my 



d 3 k 
(27r)a 



tanh K" + ~ & + tanh ^ = ^ - tanh ^ = ^ 
2 2 2 

tanh ^±±il) _ tanh ^ + "> + tanh ^t±A 



[n{uj + — /i) + n{uj- — /i) — n(& — /x) — n(w+ + /j) 
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-n(u- + /i) +n(k + n)] 



(46) 



which agrees with that of the two-loop approximately self-consistent ^-derivable HTL re- 
summation of Blaizot et al |2|, \3\. 

Now, the pressure is obtained by integrating the first line of (T4"6]) w.r.t. /i as 



V QP {T,fi) = 2N f N c T 



d 3 k 



In (1 + 



-/3(cj+-/i) 



) + h I- 



-fi(w--n) 



_|_ e -p(k-/j) 



+ In (1 + e -«^+/*)) + In I p -^^ 



(47) 



+ e _ /3 (/ c+M ) 

This agrees with the form given for quasiparticle contribution by Andersen et al J4J consider- 
ing the first term 2 in f JT5|) for /i = 0. Both quasiparticles with energies u + and cj_ generate 
T dependent ultra-violate (UV) divergences in LO HTL pressure, which is an artefact of 
1-loop HTL approximation |2j-|4j. At very high T, uj± — > k and ( 1471) reduces to free case as 
obtained in ( 1371) . 

The corresponding HTL QP entropy density in LO can be obtained as 



dv QP 



dT 



C iV/ 



d 3 k 

(2^f 



ln(l + 



-p(w+-p) 



)+.n / 



-/3(w--j0 



_|_ e -P(k-li) 



+ In (1 + e-«-+">) + In (I±^^ + fc4 + «"' " "> 



g-/8(fc+/x) y e p{u) + -[i) _|_ ]_ e /8(w_-|i) _)_ J 



(48) 



e P{k-n) _|_ X e /3(w++M) _)_ 1 e /3(w_+M) + 1 e /3(fc+M) -|_ 1 

which agrees with that of the 2-loop approximately self-consistent ^-derivable HTL resum- 
mation of Blaizot et a/ |2|. 

The QNS in LO due to HTL QP can also be obtained from (T4"6"j) as 



5 



^ ^ 5 



P? P 



4iV c iV^ 



m=o 



d 3 £; 

(2^y 



3l n(u + )(l-n(w + )) + n(w_) (1 - n(w_)) 



n(fc)(l-n(fc))] 



(49) 



where the /x derivative is performed only to the explicit /i dependence. Obviously ( 1491) agrees 
exactly with that of the 2-loop approximately self-consistent ^-derivable HTL resummation 



2 We note that the expression for QP pressure in one-loop HTLpt [4| was obtained by adding and subtracting 
the free gas pressure. However, in our formalism the correct LO form comes out naturally and no addition 



and subtraction is required as in Rcf. (J|. This is because the fluctuation of the conserved density is 
appropriately taken into consideration in the present formalism. 
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of Blaizot et al [3|. The above thermodynamical quantities in LO due to HTL quasiparticles 
with excitation energies u± are similar in form to those of free case but the hard and soft 
contributions are clearly separated out and one does not need an ad hoc separating scale as 
used in Ref. J7|. 

2. Landau Damping part (LD) 



The LD part of the quark number density follows from ( 140 p and (12"T|) as 

/"<!■,„) = N c N f / **, /^(^4W + (^) + /3-(«, mtanh^""' 



(2vr) 3 y 2tt Vw 2 - A; 2 

-k 
k 



d 3 k f f 2m 2 
N r N f I t^^ / duj I - q — ) p+(w, jfc) [n(w - aO - n(w + /i)] 



f c XV/ 



(2vr)= 



w 2 — A; 2 



(50) 



One can obtain the pressure due to LD contribution by integrating (|50p w.r.t. /x as 

k 

V LD (T, /x) = N c N f T [^ [du (^_) fl + (u, k) [In (l + e -«-/0) 



(2nY 



uj 2 — k 2 



+ ln(l + 



-/3(^+M) 



) + H 



(51) 



which has UV divergence like Andersen et al 4j and can be removed using the appropriate 
prescription therein. 

The corresponding LD part of entropy density can be obtained as 

r d 3 k 

f ' J2^f 



S LD {T^) = N c N f I ^ fdu;(J^)^(u,k) | In (1 + e-««-*>) 



+ ln (1 + e-K" + ri) + ^ - ^ + P^ + ti 
-t- yx -r e ; -r e/3(a; _ M ) + j t e ^ (w+/l) + 1 

Also the LD part of the QNS becomes 



X/ D (T)=^[p^(T,/i 



/x=0 



r d 3 k r ( 2m 2 
-it 

x /3+(o;, fc) n{u) (1 — n(a>)) , 



(52) 



(53) 



where the /x derivative is again performed only to the explicit /x dependence. It is also 
to be noted that the LD contribution is of the order of m* The LD contribution can 
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not be compared with that of the 2-loop approximately self-consistent ^-derivable HTL 
resummation of Blaizot et al 3[ as it does not have any closed form for the final expression. 
The numerical values of both the QNS agree very well. 

It is clearly evident that the various LO thermodynamic quantities in 1-loop HTLpt can 
be obtained within this modified formalism at ease instead of pushing the calculation to a 
more involved approaches. Below we demonstrate the correct inclusion of the perturbative 
content of the order g 2 to the QNS in HTLpt in a strict perturbative sense by comparing 
with the usual perturbation theory. 

V. QNS IN PERTURBATIVE LEADING ORDER (g 2 ) 



In conventional perturbation theory, for massless QCD the QNS has been calculated 16, 



18| upto order g 4 log(l/g) at /i = as 



We note that for all temperatures of relevance the series decreases with temperature and 
approaches the ideal gas value from the above, which is due to the convergence problem of 
the conventional perturbation series. 

Nevertheless, the perturbative LO, g 2 , contribution is also contained in HTL approxima- 
tion [l|, |l4j through the TV-point HTL functions. In the left panel of Fig. H] we display the 
LO HTL QNS and LO perturbative QNS scaled with free one as a function of m q /T, the 
effective strong coupling. In the weak coupling limit both approach unity whereas the HTL 
case has a little slower deviation from the ideal gas value than the LO of the conventional 
perturbative one. The latter could be termed as an improvement over the conventional 
perturbative results. The resu.ts are in very good agreement with that of Ref. & 

Next we consider a ratio |3[ as 

R= Xhtl ~ Xf , (55) 

X P (g 2 ) - Xf 

which measures the deviation of interaction of xmi from that of pQCD to order g 2 . In the 
right pane, of Fig. ffl we display this ratio as a funetion of m,/r, which approaches unity 

in the weak coupling limit indicating the correct inclusion [3[ of order g 2 in our approach 
in a strictly perturbative sense. This comes from the uj + branch [l5j of the HTL dispersion 
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N 





FIG. 4: (Color online) Left panel:The ratio of 2-flavour HTL to free quark QNS and also that of 
LO perturbative one as a function of m q /T. Right panel: The interaction measure R as a function 
of rriq/T. 

relations, u: + (k) « k + m 2 Ik, at hard momentum scale, ie, fc ~ T. With this one can now 
trivially show by expanding the QP contribution from (14911 in Sect. IV that becomes 



X? P = Xf(l-9 2 /^ 2 + •••), (56) 

which agrees with that in (j54l) . The LD contribution is of the order of mt. 

We note that the HTL resummation technique provides a consistent perturbative expan- 
sion for gauge theories at finite temperature and/or density. As discussed going beyond the 
lowest order bare perturbation theory for quark number density, we use the HTL resummed 
propagator and quark-gluon vertices in Fig. [TJ The resummed HTL quark propagators cor- 
respond to static external quarks (valence quark). In 1-loop HTLpt (viz., Fig. [T|) there is 
no dynamical quark (no quark loop) and in this sense 1-loop HTlpt is comparable with 
the quenched approximation of lattice QCD 191 ] . The inclusion of dynamical quark loops 
requires one to consider the higher-order diagrams within HTLpt in which HTL resummed 
gluon propagators (containing quark loops) will show up. This could be taken care through 
(ITTj) as it contains the covariant derivative with gauge coupling and the calculation is in 
progress. 
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VI. CONCLUSION 

In the literature the HTL resummation have been used through various approaches to 
calculate the thermodynamic quantities and also the response of the system, viz., the QNS to 
an external perturbation, i.e., the quark chemical potential. This led to different results in 
LO indicating the sensitivity of the methods. In this paper we revisited the thermodynamic 
quantities and in particular the QNS in LO within HTLpt to arrive at similar results within 
the various HTL approaches. For this purpose we modified the existing HTLpt 4j at the 
first derivative level of pressure by incorporating an infinitesimal variation to an external 
source, viz., the quark chemical potential that disturbs the system only slightly. We show 
that the various thermodynamic quantities and the QNS in LO order agree with those of 
the two-loop approximately self-consistent ^-derivable HTL resummation approach [2|, y| 
existing in the literature. Our calculation also shows that the soft and hard momenta get 
separated out naturally and one does not require any ad hoc separating scale as in [7|. 
All the thermodynamic quantities turned out to be dependent on the chemical potential 
automatically due to the method employed. We also discussed that our formalism can also 
be extended for higher-order calculations. 
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